Width distribution of contact lines on a disordered substrate 
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We have studied the roughness of a contact line of a liquid meniscus on a disordered substrate 
by measuring its width distribution. The comparison between the measured width distribution and 
the width distribution calculated in previous works, extended here to the case of open boundary 
conditions, confirms that the Joanny-de Gennes model is not sufficient to describe the dynamics of 
contact lines at the depinning threshold. This conclusion is in agreement with recent measurements 
which determine the roughness exponent by extrapolation to large system sizes. 
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The physics of elastic interfaces in random media is in- 
volved in a vast class of problems, such as domain walls 
in ferromagnetic or ferroelectric Q systems, and prop- 
agation of cracks in solids H . A notorious example of an 
elastic interface is provided by the contact line of a liquid 
meniscus on a disordered substrate 0, H H ■ 

In the past years, much effort has been devoted to shed 
light on both the equilibrium properties and the dynam- 
ics of this system [g, lJ llOj - which is characterized by 
long-range interactions. Very recently, experiments with 
water or Helium moving on a substrate characterized by 
a well controlled disorder 0, , have explored the de- 
pinning threshold. In this regime [l2l Il3| the contact 
line, driven by an external force, moves very slowly. The 
study of the roughness of the interface, and in particular 
its scaling behavior, turns out to be a fundamental tool 
to test our understanding of the physics of systems in 
which the elasticity and the disorder compete in deter- 
mining the shape of the interface. 

In Fig. we display an experimental sample: a glass 
plate with Chromium impurities (clear dots), partially 
covered by a water meniscus (dark region) . The so-called 
contact line is separating the wet and dry regions. The 
contact line is defined by an internal coordinate x and 
by the height h along the motion direction. We ob- 
serve that a single- valued height function h(x) is suffi- 
cient to characterize the shape of the contact line. A 
sample of size L — 500 ^m extracted from the image is 
also shown in Fig. ^ to analyze its geometric properties 
we are interested in the deviations from the mean height 
u(x) — h(x) — {h), with {...} = 1/L dx . . .. The mean 
square width is defined as w 2 = (u 2 (x)). The roughness 
exponent £ is introduced by considering an ensemble of 
lines of size L, where L must be larger than the lengths 
which characterize the disorder. Averaging over this en- 
semble, we obtain w 2 oc 




FIG. 1: Bottom: image of a water meniscus receding on a 
disordered substrate. The wet region is darker. Chromium 
defects appear as clear dots on the dry region. Top: sample of 
size L = 500 /im extracted from the digitalized contact line. 
u(x) is the deviation from the mean height averaged over the 
sample. 



The theoretical evaluation of the exponent £ has been 
up to now based on the assumption that the motion of the 
line at the threshold is quasi-static Eli this assumption 
was shown to be valid for most viscous fluids 11]. This 
means that the equation of motion of the height h(x) can 
be derived from an energy function, which incorporates 
the potential energy due to the driving force / and the 
disorder potential rj(x,h), as well as an elastic energy. 
According to this hypothesis, the equation of motion for 
h{x) at zero temperature is: 



d_ 

di 



h{x) = / + r)(x, h) — k dx\ 



h(x) — h{x{) 
(x — X\) 2 



(1) 
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The last term in Eq. accounts for the long-range elas- 
tic force calculated by Joanny and de Gennes |j]. At 
the equilibrium, independent approaches within Eq.^lcd 
consistently to the value ^ = 1/3 (jjlaLLJ]- At the depin- 
ning threshold, the determination of the roughness ex- 
ponent stimulated a large debate II El El 13 : finally, 



2 



extended renormalization group calculations up to the 
two- loop order proved that £ is larger than 1/3 |l7|. This 
finding was confirmed by a numerical study by means of 
an exact algorithm, which is able to detect directly the 
blocked interface at the depinning threshold: the precise 
resulting value is C = 0.388 ± 0.002 [l8|. 

In spite of the large amount of theoretical work devoted 
to the subject, experiments are very few (see and 
references therein). The main difficulty consists in con- 
trolling the disorder; in Ref. 0] this difficulty has been 
overcome by using photolithographic techniques. Small 
squares of Chromium (size: 10 x 10 /im 2 ) were deposited 
randomly on a glass plate, such that the 22% of the sur- 
face is covered. This procedure generates a disorder cor- 
related on a scale £ ~ 10 jtim, sufficiently large to prevent 
thermal fluctuations from playing any role. The corre- 
lation length £ is more than two orders of magnitude 
below the capillary length L c (~ 2.5mm in this system) 
where gravity begins to limit the fluctuations. When the 
experiment is carried out, the glass plate is withdrawn 
very slowly from the liquid bath at a fixed velocity rang- 
ing between 0.2 and 20 fixa/s. The liquid is pure water 
or an aqueous solution of glycerol with a viscosity up to 
20 times that of water. One observes that the shape of 
the contact line is independent of the velocity. This is a 
clear signature of the depinning limit [ll|. From these 
measurements it has been obtained, for 2£ < L < L c /2, 
C = 0.51 ± 0.03, in disagreement with all theoretical pre- 
dictions [Tij ]. 

The discrepancy between the theoretical and the mea- 
sured roughness exponent suggest that a richer model, 
more complicated than the one described by Eq. ^ is 
needed to account for the critical behavior of the contact 
line at the depinning. However, as the range of accessible 
scales is less than two orders of magnitude, the rigorous 
determination of the exponent £ is a very delicate experi- 
mental task. Thus, to reach more convincing conclusions, 
it is desirable to compare some other universal quanti- 
ties characterizing the depinning threshold. In particu- 
lar, special care must be devoted to evaluate finite size 
effects, which could affect drastically the final results. To 
this purpose, we study here the complete width distribu- 
tion P(w 2 ), and we compare our theoretical predictions, 
based on the solution of Eq. ^ to experimental data. 

The theoretical studies performed on systems without 
disorder such as stochastic models [13, El EH or mag- 
netic systems have shown that it is possible to express 
P(w 2 ) in a universal scaling form: 
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P(w 2 ) = =(j)(x = w 2 /w 2 ), 
w 2 



(2) 



where x is the renormalized width w 2 /w 2 . In Eq. 
the size dependence appears only through the average 
w 2 : the non-trivial scaling function <p(x) is universal and 
characterizes a full class of systems. Recently |2J], the 
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FIG. 2: Experimental scaling function (f> exp (x): study of fi- 
nite size effects. Open symbols: Aqueous solution of glyc- 
erol, v = 1 fim/s (square: L = 124 pm, circle: L = 186 /im, 
triangle: L — 372 (jm). Full symbols: water (circle: L — 
366 fim,v = 20/^m/s, square: L — 500/im,u = 1/im/s). We 
observe that for L > 300 /im the finite size effects are negli- 
gible and the scatter of the data is mostly due to the finite 
width of histogram bins. In this way we have access to the 
large scale form of <^ x P (2:). 



same scaling properties have been proved for elastic inter- 
faces in random media at the depinning threshold. The 
width distribution of these systems is, for all intents and 
purposes, given by a generalized Gaussian approximation 
of independent modes, which decay with a characteristic 
propagator G(q) ~ q d+2 ^ . Within this approximation, 
for interfaces of a fixed internal dimension and for the 
required boundary conditions, <j)(x) depends only on the 
roughness exponent £. 

In order to determine the experimental width distribu- 
tion we have employed the same experimental setup as 
described in Ref. [ll|; the contact line is imaged with 
a progressive scan CCD camera equipped with a mi- 
croscope. Two different magnifications are used, corre- 
sponding to a pixel of size Ax = 6.10 and Ax = 2.14 /xm. 
After the analysis of the images we obtain a digitized line 
of 760 pixels. Thanks to the good contrast and signal- 
to-noise ratio of the CCD, the final resolution is still one 
pixel. The experimental lines are cut into samples of size 
L = nAx, where n is the number of pixels in the sample. 
From each line we extract an ensemble of (760 — n) sam- 
ples, whose width is simply given by w 2 — 1/^X^Li u f- 
These operations are repeated over all the configurations 
detected by the camera. Clearly the values of w 2 are 
not independent ^SlJ, but this procedure ensures that no 
information is lost. 

The histogram derived by the data gives access to the 



universal function <j) Bxp (x) which is plotted in Fig. [2] for 
various sizes L and for various experimental conditions. 
For samples of size in the range 300 /im < L < 800 /im, 
the obtained function 4> exp (x) is independent of the size, 
the viscosity and the velocity of the receding meniscus. In 
smaller samples (100 fira < L < 300 /im), the distribution 
is sensitive to finite size effect related to the details of the 
way the disorder is created. Samples bigger than L = 
800 /*m cannot be treated due to the significant statistical 
noise. 

At this stage, the universal function <fi(x) which we 
have extracted from experiments can be compared to the 
theoretical evaluation within the Gaussian approxima- 
tion. Previous theoretical works EJ El Hi deal 
with samples where the boundary conditions were peri- 
odic, which is obviously not the case of our experiment 
(see Fig. 0. In order to calculate 4>{x) in the Gaussian 
approximation for open boundary conditions, we general- 
ize the discussion of Ref. l24l . The main difference, as it 
is briefly discussed in R.efs. |2.4l2fij . lays in the Fourier de- 
composition of the path u(x). In the case of open bound- 
ary conditions, the general path of size L takes the form: 



i(x) — a n cos(q n x), 



(3) 



where q n — ^ . The probability associated with this 
path is: 



(4) 



where AT is the normalization factor and, for large L, the 
exact disorder-averaged 2-point function G n takes the 
form G n — * C/n 1+2 ^. The expression for P(w 2 ) follows 
from the generating functions of the moments: 



W(z) = / dw z P{w z )e 



(5) 



Similarly to the case of periodic boundary conditions 
uM. we write 



w(z) = n 



i 



1/2 



The function <fi(x) is given by W(z) through an inverse- 
Laplace transform: 
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with A — i J2n°=i ^ an< ^ a = 1 + 2C. This complex 
integral can be written as a sum of all the tadpole con- 
tributions: 



<f>(x 



)=E^-/ ( /--n 

„=0 J ai(n) m=1 



m a A 



z + m, a A 



1/2 



, (7) 
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FIG. 3: Top: we compare the scaling function <^ox P (a;) ob- 
tained from the experimental data (squares: water L — 
500 /im, v = 1 /im/s) with the curves calculated by means of 
Eq.Q We notice the good agreement with the curve obtained 
for £ = 0.505 (continuum line). Bottom: we compare the 
scaling function 4>num(x) obtained from the numerical data 
(circles: L = 32, 10 6 independent samples) with the curves 
calculated by means of Eq. [7] We notice a perfect agreement 
with the curve obtained for £ = 0.39 (dashed line). 



with ai = -A(2n + 2) a and a 2 = -A(2n + 1) Q . The 
product in Eq. [7| converges slowly, and several thousand 
terms need to be computed. Once this is done, however, 
only a few terms in the external sum of Eq.0are sufficient 
to obtain <p(x) with high precision. 

In Ref. |3 it is discussed how the width distribution 
of an interface moving in a random medium, character- 
ized by short range elastic interactions, is described in 
an excellent way by the distribution obtained within the 
Gaussian approximation. For sake of completeness, we 
show the scaling function ^ n um(a:) computed from the 
numerical study of Eq. \T\&t the depinning threshold. We 
use the algorithm from Ref. 01 an d we calculate critical 
lines of size L — 256 with periodic boundary conditions. 
4>rmm{x) is obtained with the same procedure employed 
to deal with experimental data, by using samples of size 
L = 32. In this way, we find the scaling function corre- 
sponding to the case of open boundary conditions. 
In Fig. |31 we summarize our results. As expected from 

n derived from the numerical 



Ref. |2J], the function 
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study of Eq. ^ is in perfect agreement with the function 
obtained in the Gaussian approximation with £ = 0.39. 
The shape of 4> cxp is also in good agreement with the 
calculated function, but it is clearly shifted with respect 
to 4> nU m and best approximated with £ = 0.505. 

To our knowledge, this kind of analysis is applied for 
the first time to experimental data. It offers the possi- 
bility to evaluate the roughness exponent from a unique 
value of L and it can be useful if the accessible range of 
scaling is small. Moreover we have detected finite size 
effects that are invisible on the simple representation of 
w 2 (L). Due to the reduction of range of scaling, the anal- 
ysis of w 2 (L) yields ( = 0.52 ± 0.04, slightly larger than 
the previous determination [Til ]. 

Our study confirms that a real contact line cannot be 
described by the Joanny-de Gennes model. The same 
type of equation of motion has been proposed to describe 
the propagation of crack fronts la. In this case the ex- 
periments yield ( ~ 0.5 — 0.6 |27ll28|. Dynamical mecha- 
nism have been introduced recently |29j to account for the 
anomalous value of £ for crack fronts. On the other hand, 
for the contact line, a precise analysis of the motion justi- 
fies the quasi-static hypothesis ^lj an d no such dynamical 
mechanism has to be considered. The study of the width 
distribution of crack fronts could help to understand if 
the two systems belong to the same universality class 
and to unravel the origin of the discrepancy between the 
theoretical and the experimental values of £■ In fact, in 
both cases, the derivation of the long-range elastic term 
in Eq.Q]is obtained by a development to first order in the 
deformation. This truncation avoids the existence in the 
equation of motion of non-harmonic corrections which, 
as it has been shown for the short-range elastic interac- 
tions [33, can drastically change the critical behavior of 
the interface. Moreover the presence of non-linear terms 
has already been introduced in Ref. l Q\ to understand 
the dynamics of a contact line in the high velocity limit. 
It would be very interesting to study the effect of these 
terms on the critical properties at the depinning transi- 
tion. Unfortunately up to now the numerical computa- 
tion of this new equation of motion is impossible because 
of the presence of non-convex terms in the elastic energy. 

We thank C. Guthmann, P. Le Doussal, J. Vannimenus 
and J. K. Wiese for useful discussions. 
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